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We determine the decoupling surfaces for the hydrodynamic description of heavy ion collisions at 
RHIC and LHC by comparing the local hydrodynamic expansion rate with the microscopic pion- 
pion scattering rate. The pion pt spectra for nuclear collisions at RHIC and LHC are computed 
by applying the Cooper-Frye procedure on the dynamical-decoupling surfaces, and compared with 
those obtained from the constant-temperature freeze-out surfaces. Comparison with RHIC data 
shows that the system indeed decouples when the expansion rate becomes comparable with the pion 
scattering rate. The dynamical decoupling based on the rates comparison also suggests that the 
effective decoupling temperature in central heavy ion collisions remains practically unchanged from 
RHIC to LHC. 



INTRODUCTION 

Non-dissipative hydrodynamical description has been 
quite successful in explaining the data on low-p^ trans- 
verse momentum spectra and elliptic flow in nucleus- 
nucleus collisions at RHIC P, i, [1 i, i, S, 0, i- First 
results from viscous hydrodynamics suggest, that the vis- 
cosity [13, U, [I2II of quark-gluon plasma should be 
small to accommodate the data, in particular those on 
elliptic flow. Apart from the viscous effects, the main un- 
certainties in the hydrodynamic modelling of heavy ion 
collisions lie in the computation of the initial conditions, 
the choice of equation of state (EoS) for the QCD mat- 
ter and, the subject of this paper, the treatment of the 
freeze-out from a local thermal equilibrium state to free 
hadrons. The standard way of handlin g th e freeze-out is 
the Cooper-Frye decoupling procedure p3| . Also hadron 
cascade models have been successfully implemented to 
describe late hadronic dynamics [IJ, |l5|, but typically 
the initial hadron distributions in these simulations are 
also calculated using the Cooper-Frye method. 

Traditionally, the transverse momentum distributions 
of hadrons at freeze-out are calculated as emissions from 
a constant-temperature hypersurface, with the value of 
the decoupling temperature Tdoc fixed from data. Since 
the constant-Tdcc Cooper-Frye decoupling procedure is 
not based on scattering rates, it cannot predict the freeze- 
out temperature or conditions for the LHC Pb+Pb col- 
lisions. To obtain predictions for hadron pT spectra at 
the LHC, one often assumes - like we have done in 0] - 
that Tdcc at RHIC and LHC are equal. 

In this work, we shall study this uncertainty by de- 
termining the decoupling surface by using a dynamical 
freeze-out condition instead of a constant Tdcc- That 
is, we compare the local hydrodynamical expansion rate 
with the local scattering rate of pions, and define the 
Cooper-Frye freeze-out surface to be where the expan- 
sion rate exceeds the scattering rate, possibly multiplied 
by a constant. To our knowledge, this idea was first dis- 
cussed in Refs. 3 12|- Impact parameter dependence 



of freeze-out conditions using the same idea is studied in 
Ref. [18(1 . The ratio of these two rates can be fixed on the 
basis of the measured pion px spectra at RHIC. Because 
this ratio now depends only on the local properties of the 
system (temperature and flow) it can well be expected to 
remain unchanged in the extrapolation to higher collision 
energies. We are thus able to predict the freeze-out sur- 
face for the bulk of QCD matter in Pb-|-Pb collisions at 
the LHC based on the dynamics of the system. 

We calculate pion pT spectra using freeze-out surfaces 
from the dynamical decoupling condition, and compare 
the results with the constant-T decoupling at RHIC and 
LHC. Interestingly, two main results come out of this 
exercise: First, comparison with the pion pT spectra at 
RHIC nicely confirms the expectation that the expansion 
and scattering rates are of the same order at freeze-out. 
Second, the local dynamical freeze-out condition leads to 
practically the same effective decoupling temperature at 
LHC as at RHIC which is non-trivial since the evolution 
of flow is quite different at the two collision energies. 



HYDRODYNAMIC SETUP AND RATES 

A detailed account of the hydrodynamic setup we use 
can be found in Q- We consider boost-invariant per- 
fect fluid hydrodynamics with transverse expansion for 
azimuthally symmetric systems, i.e. study (nearly) cen- 
tral collisions only. The EoS is given by the bag model, 
where the quark-gluon plasma (QGP) is an ideal gas of 
massless quarks, antiquarks and gluons, and the hadron 
resonance gas (HRG) Il9l consists of all hadronic states 
up to a mass 2 GeV [20|. Critical temperature is fixed 
to Tc = 165 MeV at ^jlb =0. Full thermalization, both 
kinetic and chemical, until decoupling is assumed, after 
which we also include the 2- and 3-body resonance decays 
of unstable hadronic states. 

The initial conditions for the hydrodynamic evolution 
- the total energy and net-baryon number deposited into 
a rapidity unit by an initial time - for RHIC and LHC 
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A + A collisions, are obtained from the pQCD + (final 
state) saturation model j2l||as explained in detail in 
In our previous studies [j, |7[ the initial energy density 
was taken to be proportional to the density of nucleon 
binary collisions (BC). To probe the uncertainty due to 
the initial transverse geometry, we study also an initial 
state where the energy density scales with the density of 
wounded nucleons (WN). In the latter case, the initial 
entropy is kept the same as in the BC case. 

A measure of the hydrodynamical expansion rate is 
given by the expansion scalar 
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where m'' — 7(1, w^, w^, — 7r(cosh77, w^, 0, sinhiy) is 
the flow four-velocity of the matter, rj is the coordinate 
space rapidity and r the longitudinal proper time. This 
can be interpreted as the rate of relative volume change 
V/V in the local rest frame [l6[, with V = dV/dr' , where 
t' is the fluid element's proper time. An easy way to see 
this is to express entropy conservation in the local rest 
frame as 

s = -s0, (2) 
where s is the entropy density. From s = S/V it follows 

(3) 



V 



where the last step follows from S = constant. The ex- 
pansion rate 6 can be easily calculated once the hydro- 
dynamic solution is known. 

As long as the microscopic collision rate per particle, 
r, is large relative to the expansion rate, T ^ 6, the sys- 
tem can be assumed to be in local thermal equilibrium. 
Since we are here interested in the bulk matter proper- 
ties only, and since without studying a separate chemical 
decoupling of heavier particles we cannot reproduce their 
Pt spectra in detail [n], we consider only pion-pion col- 
lisions in the scattering rate. We thus approximate the 
local collision rate per particle as 



(4) 



where the density of (massless) pions is nTj{T) — 2)^^T^ . 
For the average pion-pion cross-section at a given tem- 
perature r, we use the following parametrization, 



C^TTTriT) = (To 



(T- 



(A/2)2 



(5) 



where ctq = 0.60 fm^, C = 0.78 fm^, To = 105 MeV and 
A = 170 MeV, which reproduces the results of Ref. [23| 
in the temperature range T < 165 MeV of the HRG 
phase. Note in particular, how due to the rho resonance 
peak the cross section decreases with decreasing T, in- 
stead of a naive dimensional T-scaling cr^jr ~ T~^. Note 



also that we do not need to consider partonic scatter- 
ing rates at all. Since we are interested in the decou- 
pling of final hadrons, we can limit the rate study to the 
HRG phase only. In the Bjorkcn picture [2^ of heavy ion 
collisions, which is adopted here, the hydrodynamic ex- 
pansion starts with very strong longitudinal expansion. 
In the boost-invariant case without transverse expansion, 
the fluid velocity, = (cosh 77, 0, 0, sinhjy), depends only 
on the coordinate space rapidity 77, and the (longitudinal) 
expansion rate is given by the inverse of the longitudinal 
proper time t, 9 — 1/t. During the expansion, pres- 
sure gradients develop a strong transverse flow, which 
leads to an increasing transverse expansion rate. The 
decreasing longitudinal and the increasing transverse ex- 
pansion rates lead to a system where the local expansion 
rate is high at the spatial boundary of the system, where 
transverse gradient is strong, while near the center of the 
fireball the expansion rate is relatively weak. 
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FIG. 1: Contours of local expansion rate 6, phase boundaries 
and two isotherms in 5 % most central RHIC 200 GeV Au-|-Au 
collision for BC initial state. Reading from the top of the 
figure, 9 = 0.2, 0.3, . . . 1.8. Also two flowlines, starting from 
r = 3 fm and r — 6 fm, are shown. 

In Fig. 1 we show the contours of constant 6 for a nearly 
central RHIC Au-I-Au collision at ^^snn = 200 GeV 
together with the phase boundaries and two isotherms 
T = 150 MeV and T = 120 MeV. Of the two isotherms 
shown, the one with T = 150 MeV corresponds to the 
freeze-out surface, which leads to an agreement with the 
RHIC data on pion px spectra, as discussed in 7] . Also 
two flowlines starting from transverse locations r = 3 fm 
and r = 6 fm are drawn in the figure. 

Figure 2 shows the computed expansion and scattering 
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FIG. 2: Expansion and scattering rates along the flowlines 
starting from r = 0, 3 and 6 fm for RHIC 200 GeV Au+Au 
collision. Flowlines are shown in Fig. 1. Also the expansion 
rate 9 = 1/t for a system with longitudinal boost invariant 
expansion only is shown. 

rates at r = fm and along the two flowlines in Fig. 1. 
For comparison with the case of no transverse expansion, 
also the expansion rate 6 = 1/t, due to the longitudinal 
boost invariant flow only, is shown. As seen in Fig. 2, the 
system starts its evolution with longitudinal but without 
transverse expansion and 9 drops like 1/t until the in- 



creasing transverse expansion starts to contribute to the 
rate, slowing down its decrease. Interestingly, this can 
even develop a local maximum of 9 near the boundary 
of QGP and mixed phase. In Fig. 1 the rarefaction wave 
from the system boundary can be seen as a bending of 
the 9 contours towards the center of the fireball in the 
QGP phase. The same phenomenon can be seen in Fig. 2 
when the expansion rate starts to exceed the longitudinal 
expansion rate. 

Because pressure gradients are stronger near the 
boundary of the system, the transverse expansion rate 
grows faster at larger values of transverse radius r. This 
is seen clearly in Fig. 1, where in the QGP phase the 
expansion rate grows as a function of r at constant r. In 
the mixed phase pressure gradients vanish, which stops 
the increase of the expansion rate. In contrast to the 
behaviour in the QGP phase, in the mixed phase the ex- 
pansion rate decreases as a function of r at constant r. 
This is because the matter has spent less time in the QGP 
phase at larger values of r, thus there has been less time 
for transverse flow to develop. When the system goes 
through the phase transition into HRG, the shock wave 
formed generates very strong gradients near the phase 
boundary in the HRG, and this in turn leads again to a 
strong increase of the transverse expansion as a function 
of r in the HRG phase. 

In summary, the phase transition first slows down the 
increase of the transverse expansion rate, but this is later 
counteracted by resulting strong gradients in hadron gas 
near the phase boundary. This is seen in Fig. 1, where, as 
a result of the phase transition, the 9 contours are com- 
pressed near the mixed phase and hadron gas boundary. 
This effect will be important when we construct the de- 
coupling boundaries by comparing the local scattering 
rate to the local expansion rate. 

In Fig. 2, for the flowline starting from r = 6 fm (see 
the lowest panel) , we see that 9 can actually decrease be- 
low the longitudinal value 1/t for a short while during the 
phase transition. This is due to a transverse compression 
rather than expansion in a small region near the bound- 
ary during the phase transition, caused by the QGP ac- 
celerating against the mixed phase where the pressure 
gradients are zero and matter is not accelerating. This is 
how compressional shock wave is mechanically formed be- 
tween the QGP and the HRG phases. Same phenomenon 
can be seen in Fig. 1, where expansion rate has a strong 
local minimum in the mixed phase near the edge of the 
fireball. 

Because the scattering rate is a function of tempera- 
ture alone, the constant-temperature contours represent 
also the contours of the scattering rate. It is clear from 
Fig. 1 that the contours of expansion and reaction rates 
do not coincide. Both the scattering and expansion rates 
decrease in general with time, but the decrease of the 
scattering rate is faster - thanks to the resonance peak 
in the pion scattering cross section. At sufficiently low 
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temperatures, this eventually leads to the situation where 
the scattering rate drops below the expansion rate, as il- 
lustrated in the two top panels of Fig. 2, and this is when 
the system is expected to decouple. Thus, the decoupling 
surface is given by the crossing points of the expansion 
and scattering rate lines like in the two top panels of 
Fig. 2. 

The bottom panel of Fig. 2 shows the situation near 
the edge of the system, where the scattering rate in the 
HRG falls below the expansion rate already at the be- 
ginning of the HRG phase. In this case, we assume the 
decoupling to take place at the boundary of the mixed 
phase and the HRG. 

Figure 3 shows the phase boundaries and constant tem- 
perature contours as in Fig. 1, and in addition also the 
decoupling boundaries from a dynamical condition 



Au + Au 0-5% 



e = cF, 



(6) 



with different values of the constant c. (The crossing 
points in Fig. 2 correspond to c = 1.) In the hydro- 
dynamic approach with phase transition we cannot de- 
couple the matter before it has completely hadronized. 
Thus the decoupling curve can be only in HRG or on the 
boundary of mixed phase and HRG. Dynamics of the sys- 
tem with phase transition is such that it always forces de- 
coupling curves near the mixed phase-HRG phase bound- 
ary, as seen in Fig. 1, where contours of expansion rate are 
compressed near the boundary. We only need a slight ad- 
justment to force the freeze-out boundaries to be strictly 
outside the mixed phase. It turns out that for the values 
of c considered here, the rate of transverse expansion on 
the edge of the matter is so large that the decoupling 
curve always first follows the boundary of the mixed and 
HRG phases. For the smallest value, c = 0.5, the decou- 
pling curve (see Fig. 3) is close to the mentioned phase 
boundary at all radial distances, but for larger values of 
c the curves separate when the phase boundary curve 
turns sharply towards the center of the matter. Near the 
center, where the expansion rate is smaller than at the 
edge since the contribution from the transverse expan- 
sion is small, the freeze-out time can become large and 
the freeze-out temperature small. Thus along a dynamic 
decoupling curve the temperature starts from = 165 
MeV at the edge of the matter and drops, for c>l to 
clearly smaller values of the order of 100 to 120 MeV 
as the curve reaches the center. As seen in Fig. 3, the 
dynamical decoupling curves typically cross the constant- 
temperature curves. 

For large values of c the rather big differences in tem- 
perature along the decoupling curve can cause the follow- 
ing problem. Without a chemical freeze-out, at or soon 
after the formation of hadron phase, the density of heav- 
ier particles can become much smaller on the flat, space- 
like part of the decoupling curve than on the time-like 
part at the edge. Since in the Cooper-Frye prescription 
the contribution to the spectrum from the time-like part 
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FIG. 3: Decoupling boundaries from a dynamical condition 
d = corresponding to 5 % most central RHIC 200 GeV 
Au+Au collisions. Calculation is initialized with BC profile. 
Also phase boundaries and constant temperature decoupling 
boundaries are shown. 



can become negative at small pt's, it can happen that too 
few heavy particles are available on the low temperature 
space-like part to exceed the negative contributions, and 
the spectrum may remain negative close to = 0. From 
our point of view there is no reason to go to large values 
of c since c '--^ 1 is the natural assumption. Also, we are 
considering only pions in this work, and the described 
effect is much weaker for pions. 

Finally, we note that so far all the figures shown have 
been for the BC initial state, but we have checked that 
similar results also hold for the WN profile. 



RHIC RESULTS 

In Fig. 4 we show our previous results with BC ini- 
tial profile 0] for 5 % most central Au+Au collisions at 
y/sJfN = 200 GeV, compared with the RHIC data from 
BRAHMS ;24!j, PHENIX 25] and STAR ^] collabora- 
tions. With the WN profile for the initial energy and 
net-baryon number densities, the pressure gradients are 
smaller and the transverse flow builds up more slowly. 
Thus a lower Tdcc is needed with the WN initial pro- 
file in order to reproduce the RHIC data. The figure 
demonstrates that both the BC initial conditions with 
freeze-out at Tdcc = 150 MeV and the WN initial con- 
ditions with freeze-out at Tdcc = 140 MeV are in good 
agreement with the data and with each other. We shall 
use these constant- Tdcc results as a comparison baseline 
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FIG. 4: Transverse momentum spectra of positive pions 
at y ~ in 5 % most central RHIC 200 GeV Au+Au 
collisions, computed from constant-temperature decoupling 
boundaries for different initial states. The data shown are 
from BRAHMS [H], PHENIX [H and STAR [H] collabora- 
tions. 




12 3 4 

Pt [GeV] 



FIG. 5: Transverse momentum spectra of positive pions at 
2/ in 5 % most central RHIC 200 GeV Au-|-Au collisions 
from the dynamical decoupling condition for different values 
of the constant c, compared with the constant-temperature 
spectra. The calculation is made for the BC initial conditions. 

in what follows. Note how the freeze-out at Tdcc = 120 
MeV clearly overshoots the data, for both initial condi- 
tions. 

Figure 5 shows the p^-spectra of positive pions, ob- 
tained with the BC initialization and using the dynani- 
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FIG. 6: As Fig. 5 but for WN initial conditions. 



ical decoupHng Eq. The BC-initialized Tdcc — 150 

and 120 MeV results from Fig. 4 are shown for compar- 
ison. Figure 6 shows the corresponding results with the 
WN initialization. We notice that for the BC initial con- 
dition the c values in the range 1 — 1.5 give a reasonable 
description of the data, the spectra being all close to our 
previous calculation for Tdcc = 150 MeV. A closer inspec- 
tion indicates that c ~ 1 gives perhaps the best descrip- 
tion. The fact that the constant-Tdoc leads to very similar 
spectra as the rate condition with c ^ 1, which we con- 
sider to be physically better motivated, shows that also 
the constant-Tdcc on the average correctly captures the 
effects of the rates competition. As indicated by Fig. 6, 
the use of WN profile supports slightly larger values of 
c, with c ~ 1 — 1.5 giving the best description of the 
data at pt^2 GeV, i.e. in the most relevant region for 
hydrodynamical models. On the basis of RHIC data we 
can thus conclude that within our dynamical freeze-out 
condition, a good description of the data is obtained with 
c~ 1 - 1.5. 



LHC RESULTS 

In constructing the hypersurface for the Cooper-Frye 
decoupling based on the dynamical freeze-out condition, 
we assume that freeze-out dynamics is governed by local 
conditions only. Thus we expect that the constant c ob- 
tained above is independent of the global properties of 
the system, such as the system size or the centrality of 
the collision. In particular, we expect c to be indepen- 
dent of the collision energy. Therefore we can use the 
value of c, which is determined from the RHIC data, to 
predict the decoupling boundary at the LHC. 
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FIG. 7: As Fig. 3 but for 5 % most central LHC 5500 GeV 
Pb+Pb collision. 

Figure 7 shows the phase boundaries and decoupling 
curves for 5 % most central Pb+Pb colhsions at the LHC 
energy ^^snn = 5.5 TeV, in the case of the BC profile for 
the initial energy density. The behaviour of the decou- 
phng boundaries, which are obtained from the dynamical 
decoupling condition, is very similar to the RHIC results 
above: during the early times of the evolution decoupling 
near the spatial boundary of the system occurs immedi- 
ately after hadronization is complete, while at the center 
of the fireball the system decouples at later times, when 
decoupling temperature is significantly lower. The de- 
coupling boundaries for the WN profile are qualitatively 
similar to the BC case, and therefore not shown. 

Transverse momentum spectra of positive pions in 5 
% most central Pb-|-Pb collisions at the LHC, calculated 
from the decoupling boundaries above are shown in Fig. 8 
for the BC and Fig. 9 for the WN initializations. Fix- 
ing the proportionality constant c between the expansion 
and scattering rates from the RHIC data in Fig. 4, our 
prediction for the positive pion px spectra is given by the 
curves with c = 1 — 1.5. Also the constant-Tdcc results 
with Tdcc — 150 and 120 MeV are shown for the BC ini- 
tiahzation in Fig. 8, and with Tdcc = 140 and 120 MeV 
for WN initialization in Fig. 9. We make the following 
interesting observations: 

(i) Once the value of c is fixed from RHIC data, take 
e.g. c = 1 with the BC initialization (Figs. 4 and 8), the 
dynamical decoupling results both at RHIC and LHC are 
reproduced with the same Tdoc = 150 MeV. This sug- 
gests that Tdcc does not change considerably from RHIC 
to LHC. The same holds for the results for WN initializa- 
tion, with effective Tdcc — 140 MeV. Upper limit c = 1.5 
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FIG. 8: As Fig. 5 but for 5 % most central LHC 5500 GeV 
Pb-|-Pb collision. 



lO-" 



10-^ 



> 10' 
O 



■7^ 10° 

Q. 
■D 



i 10" 



10"' 



10"' 



10"' 





Pb + Pb 0-5% 
s^'^ = 5500 A GeV 
WN initialization 



1 2 3 4 5 6 

Pt [GeV] 

FIG. 9: As Fig. 6 but for 5 % most central LHC 5500 GeV 
Pb+Pb collision. 

corresponds a slightly smaller effective Tdoc in both cases. 
Thus, in our hydrodynamical framework, we expect ef- 
fective Tdcc to be in the range 140 - 150 MeV. 

(ii) The LHC pion spectra obtained with c = 1 — 1.5 
or with an effective Tdcc = 150 MeV for the BC case in 
Fig. 8, or with an effective Tdcc = 140 MeV for the WN 
case in Fig. 9, are all very close to each other. This 
shows that our LHC predictions for pion spectra 0] are 
not very sensitive to the uncertainties of the initial pro- 
file, provided that the RHIC data have been exploited to 
constrain the decoupling dynamics. 
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(in) As seen in Figs. 4-6, the RHIC data outrules the 
values of c which would be much larger than 1-2. Corre- 
spondingly, as shown before [1,01 , a low effective constant 
Tdec, cannot be used at RHIC. The present study of local 
dynamic decoupling condition lends support to the pre- 
diction that if constant-Tdcc decoupling is used, low Tjoc 
values are unlikely also at the LHC. 

(iv) Intuitively, since the flow develops faster already 
in the plasma phase at the LHC than at RHIC, one might 
anticipate the decoupling at the LHC to happen at higher 
effective temperatures than at RHIC. The present study, 
especially the BC case, now shows that the effective de- 
coupling temperature in central Pb-|-Pb collisions at the 
LHC is in fact quite close to that in central Au-|-Au col- 
lisions at RHIC. 



CONCLUSIONS 

We have computed the transverse momentum spec- 
tra of pions in nearly central Au-|-Au collisions at RHIC 
and Pb+Pb collisions at the LHC within the framework 
of ideal hydrodynamics, by taking the initial conditions 
from a pQCD + saturation model. In particular, we have 



where 



applied a dynamical decoupling condition 
we compare the local expansion rate to the local pion- 
pion scattering rate. The proportionality constant be- 
tween the two rates at decoupling is extracted from the 
RHIC data on pion transverse momentum spectra. The 
same proportionality constant is then used to present pre- 
dictions of pion spectra at the LHC. 

Although the model is simple and we consider pion- 
pion collisions and pion spectra only, the fact that the 
proportionality constant obtained is of the order of unity, 
gives us some confidence that we have captured an essen- 
tial part of the physics in the breakdown of the hydro- 
dynamic stage of matter into final free hadrons. The 
dynamical decoupling condition does not necessarily in- 
dicate the end of all hadronic interactions, but for tttt 
interactions this is expected to be the case, at least in 
the sense that interactions cannot maintain local ther- 
mal equilibrium. Nucleons most likely decouple some- 
what later than pions because of the strong A resonance 
in ttN interactions 14l. Il7l|. For pions, however, this is 



only a minor effect, due to the small number of nucleons 
compared to pions. A consistent modelling of these inter- 
esting details is beyond the scope of the present paper. 

Studying both binary collision and wounded nucleon 
profiles for the initial energy densities, we have compared 
the results obtained with the dynamical decoupling crite- 
rion to those obtained with effective constant freeze-out 
temperatures. We find that the effective freeze-out tem- 
peratures for pions at the LHC will be very close to those 
at RHIC, Tdoc = 140 - 150 MeV. With the dynamical 
decoupling criterion, we also find that the hydrodynamic 
prediction of the pion pT spectra at the LHC is fairly 



insensitive to the details of initial transverse profile, once 
the total entropy is fixed and once the decoupling dy- 
namics is first constrained by the RHIC data. 

We thank P. Huovinen for discussions and the 
Academy of Finland, Project 115262, for financial sup- 
port. 
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